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Detailed experimental data for physisorption potential-energy curves of H2 on low-indexed faces 
of Cu challenge theory. Recently, density-functional theory has been developed to also account for 
nonlocal correlation effects, including van der Waals forces. We show that one functional, denoted 
vdW-DF2, gives a potential-energy curve promisingly close to the experiment-derived physisorption- 
energy curve. The comparison also gives indications for further improvements of the functionals. 



PACS numbers: 71.15.-m, 73.90.+I, 68.35.Np 

Density-functional theory (DFT) gives in principle ex- 
act descriptions of stability and structure of electron sys- 
tems, but in practice approximations have to be made to 
describe electron exchange and correlation (XC) [U [2]. 
Evaluation of XC functionals is commonly done by 
comparison with results from other accurate electron- 
structure theories or by comparing with relevant exper- 
iments, typically providing test numbers only for one or 
two measurables. This Letter illustrates the advantages 
of a third approach, which builds the experiment-theory 
calibration based on extensive experimental data, in this 
case a full physisorption potential derived from surface- 
physics measurements. 

In the physisorption regime, resonant elastic 
backscattering-diffraction experiments provide a de- 
tailed quantitative knowledge. Here, data obtained for 
H 2 and D 2 on Cu surfaces [HHS] are used as a demanding 
benchmark for the performance of adsorbate potential- 
energy curves (PECs) calculated with a nonempirical 
theory for extended systems. Density functionals that 
aspire to account for nonlocal electron-correlation ef- 
fects, including van der Waals (vdW) forces, can then be 
assessed. In particular, we study the vdW-DF method 
[THTU] and show that calculations with versions of it 
provide a promising description of the physisorption 
potential for H2 on the Cu(lll) surface and that the 
most recent one, vdW-DF2 10J, is more accurate than 
the first one [THS] and other tested functionals. 

Sparse matter is abundant. Dense matter, also abun- 
dant, is since long successfully described by DFT. The 
recent extensions of DFT functionals to regions of low 
electron density, where the ubiquitous vdW forces are 
particularly relevant, render DFT useful also for sparse 
matter. In the vdW-DF functional the vdW interactions 
and correlations are expressed in the density n(r) as a 
truly nonlocal six-dimensional integral [7HH]. Its key in- 



gredients are (i) its origin in the adiabatic connection for- 
mula |11H13| . (ii) an approximate coupling-constant in- 
tegration, (iii) the use of an approximate dielectric func- 
tion in a single-pole form, (iv) which is fully nonlocal 
and satisfies known limits, sum rules, and invariances, 
and (v) whose pole strength is determined by a sum rule 
and whose pole position is scaled to give the approxi- 
mate gradient-corrected electron-gas ground-state energy 
locally. There are no empirical or fitted parameters, just 
references to general theoretical criteria. 

Like composite molecules, adsorption systems have 
electrons in separate molecule-like regions, with exponen- 
tially decaying tails in between. Then the slowly varying 
electron gas, used in the original vdW-DF method [7J- 
Hl [H] ) might not be the most appropriate reference sys- 
tem for the gradient correction |15j . Although promis- 
ing results have been obtained for a variety of systems, 
including adsorption [THl \T7\ , there is room for improve- 
ments. For the recent vdW-DF2 functional, the gradient 
coefficient of the B88 exchange functional [TB] is used for 
the determination of the internal functional [Eq. (12) of 
Ref. [7] within the nonlocal correlation functional. This 
is based on application of the large- N asymptote [HI [20] 
on appropriate molecular systems. Using this method, 
Elliott and Burke [3T] have shown from first principles 
that the correct exchange gradient coefficient j3 for an 
isolated atom (monomer) is essentially identical to the 
B88 value, which had been previously determined empir- 
ically [18 . Thus in the internal functional, vdW-DF2 [TO] 
replaces Z a b in that equation with the value implied by 
the (3 of B88. This procedure defines the relationship 
between the kernels of vdW-DF and vdW-DF2 for the 
nonlocal correlation energy. Like vdW-DF, vdW-DF2 is 
a transferable functional based on physical principles and 
approximation and without empirical input. 

The vdW-DF method also needs to choose an over- 
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all exchange functional to obtain the exchange contribu- 
tion to the interaction energy between two monomers (for 
example). The original vdW-DF uses the revPBE [22] 
exchange functional, good at separations in typical 
vdW complexes [?H5]. The latter choice can be im- 
proved on [2"3H2"T] . Recent studies suggest that the 
PW86 exchange functional [5S] most closely duplicates 
Hartree-Fock interaction energies both for atoms [24] and 
molecules [25]. The vdW-DF2 functional jTU] employs 
the PW86R functional [2~5] . which more closely repro- 
duces the PW86 integral form at lower densities than 
those considered by the original PW86 authors. 

Evaluation of XC functionals with respect to other 
theoretical results is often done systematically, e.g., by 
benchmarking against the S22 data sets [2T)H53] . The 
S22 sets have twenty-two prototypical small molecular 
duplexes for noncovalent interactions (hydrogen-bonded, 
dispersion-dominated, and mixed) in biological molecules 
and provides PECs at a very accurate level of wave- 
function methods, in particular the CCSD(T) method. 
However, by necessity, the electron systems in such sets 
have finite size. The original vdW-DF performs well 
on the S22 dataset [2"5H33] , except for hydrogen-bonded 
duplexes (underbinding by about 15% [TO] US]). The 
vdW-DF2 functional reduces the mean absolute devia- 
tions of binding energy and equilibrium separation signif- 
icantly [TO] . Shapes of PECs away from the equilibrium 
separation are greatly improved. The long-range part 
of the vdW interaction, particularly crucial for extended 
systems, has a weaker attraction in the vdW-DF2, thus 
reducing the error to 8 meV at separations 1 A away from 
equilibrium [10] . 

Experimental information provides the ultimate basis 
for assessing functionals. The vdW-DF functional has 
been promising in applications to a variety of systems 
[16], but primarily vdW bonded ones, typically tested 
on binding-energy and/or bond-length values that hap- 
pen to be available. The vdW-DF2 functional has also 
been successfully applied to some extended systems, like 
graphene and graphite [10] , metal-organic- frameworks 
systems [34], molecular crystal systems [35], physisorp- 
tion systems [3J3 13Z] ; liquid water [35] and layered oxides 
[35] , However, those studies are of the common kind that 
focus on comparison against just a few accessible obser- 
vations. 

The key step taken by the present work is to bench- 
mark a full PEC in an extended system. Fortunately, 
for almost two decades, accurate experimental values for 
the eigenenergies of H2 and D2 molecules bound to Cu 
surfaces [3J U] have been waiting for theoretical account 
and assessment. The rich data bank covers results for 
the whole shape of the physisorption potentials. 

The H2-CU system is particularly demanding for the 
vdW-DF and vdW-DF2 functionals and alike. On one 
hand, H2 is a small molecule with a large HOMO-LUMO 
gap, far from the low-frequency polarization modes as- 



sumed in derivation of these functionals [TUTU]. On the 
other hand, Cu is a metal, which has created particular 
concerns [40] . 

Chemically inert atoms and molecules adsorb physi- 
cally on cold metal surfaces [5], with characteristic low 
desorption temperatures ranging from a few K (He) to 
tens of K (say Ar and CH4). Adsorption energies, which 
take values from a few meV to around 100 meV, may 
be determined from thermal-desorption or isosteric-heat- 
of-adsorption measurements. For light adsorbates, like 
He and H2, gas-surface-scattering experiments, involving 
resonance structure of the elastic backscattering, provide 
a more direct and elegant method, with accurate and 
detailed measurements of bound-level sequences in the 
potential well. The availability of isotopes with widely 
different masses ( 3 He, 4 He, H 2 , D 2 ) permits a unique as- 
signment of the levels and a determination of the well 
depth and ultimately a qualified test of model poten- 
tials [6]. 

The bound-level sequences of specific concern here 
were obtained using nozzle beams of para-H 2 and normal- 
D2, i.e., the beams are predominantly composed of j = 
molecules. This implies that the measured bound-state 
energies, e„ (listed in Table[TJfor H 2 and D 2 on Cu(lll)), 
refer to an isotropic distribution of the molecular ori- 
entation. For this particular ordering, all e„ values fall 
accurately on a common curve when plotted versus the 
mass-reduced level number 77 = (ri + l/2)/y / m. This im- 
plies a level assignment that is compatible with a single 
gas-surface potential for the two hydrogen isotopes [4]. 
A third-order polynomial fit to the data yields for 77 = 
a potential-well depth D = 29.5 meV. 

The experimental energy levels in the H2-Cu(lll) PEC 
(see Table|T| may be analyzed [3H5] within the traditional 
theoretical picture [UJ 02] of the interaction between in- 
ert adsorbates and metal surfaces: The PEC is then ap- 
proximated as a superposition of a long-range vdW at- 
traction, K,dw, an d a short-range Pauli repulsion, Vr, 
reflecting overlap between tails of metal Bloch electrons 
and the adparticle's closed-shell electrons [UJ . This 
results in a laterally and angularly averaged potential 
V (z) — Vu(z) + V v dw( z ), where z is the normal dis- 
tance of the H 2 bond center from the jellium edge. The 
bound-level sequences in Table [TJ can be accurately re- 
produced (< 0.3 meV) by such a physisorption poten- 
tial [1 [5] (Fig. [T]), having a well depth of 28.9 meV and 
a potential minimum located 3.50 A outside the topmost 
layer of copper ion cores. From the measured intensities 
of the first-order diffraction beams, a very small lateral 
variation of the H2-Cu(lll) potential can be deduced, 
~ 0.5 meV at the potential-well minimum. 

A direct solution of the Schrodinger equation in V (z) 
reproduces the four low-energy eigenvalues to within 3% 
of the measured ones. It is therefore consistent [33] to 
benchmark against the very accurately constructed ex- 
perimental physisorption curve in Fig. [l] 
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FIG. 1: Experimentally determined effective physisorption 
potential for H2 on Cu(lll) 3 , compared with potential- 
energy curves for H2 on Cu(lll), calculated for the atop site 
in GGA-revPBE, GGA-PBE, vdW-DF2, and vdW-DF. 



TABLE I: Sequences of bound-state energies for H2 and D2 
on Cu(lll). DFT eigenvalues are calculated with the vdW- 
DF2 potential of Fig. [I] and experimental numbers are from 
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e„ [meV] 




H 2 D 2 




DFT Exper. DFT Exper. 





-32.6 -23.9 -34.4 


1 


-21.3 -15.5 -26.0 -19.0 


2 


-12.1 -8.7 -18.7 -12.9 


3 


-5.4 -5.0 -12.4 -8.9 


4 


-7.4 -5.6 


5 


-3.5 -3.3 



Figure [T] shows our comparison of density functional 
PECs against the experimental physisorption potential 
for H2 in an orientation above an atop site on the Cu(lll) 
surface. The PECs of Fig.[l] are calculated with the vdW- 
DF [7] and vdW-DF2 functional [10] as well as with 
two generalized gradient approximations (GGA-PBE and 
GGA-revPBE). We use an efficient vdW algorithm [H] 
adapted from SIESTA's [12 vdW code within a modi- 
fied version of the plane- wave code ABINIT 46J. The 
vdW interaction is treated fully self-consistently [9] (al- 
lowing also vdW forces to relax the adsorption geome- 
try). The computational costs are the same with vdW- 
DF and vdW-DF2. Our choice of Troullier-Martins-type 
norm-conserving pseudopotentials and a high cutoff en- 
ergy (70 Ry) ensures excellent convergence; the GPAW 
code [47] , in its default mode, gives similar but less accu- 
rate results [IS] . We stress that there is neither a damp- 
ing nor saturation function in vdW-DF and vdW-DF2 
calculations. The need for an account of nonlocal corre- 
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FIG. 2: Interaction potential for H2 on Cu(lll), calculated 
self-consistently with the vdW-DF2 functional |10j in the 
bridge, hollow and atop sites. 



lations for the description of vdW forces is illustrated by 
the GGA curves giving inadequate PECs. The calculated 
well depth in vdW-DF, 53 meV, should be compared with 
the measured one, 29.5 meV g], and the one calculated 
from V (z), suitably parametrized, 28.9 meV [3] (Fig. [I]). 
We find that the vdW-DF2 PEC lies close to the exper- 
imental physisorption potential, both at the equilibrium 
position and at separations further away from the sur- 
face. 

Calculated PECs of H2 above bridge, atop, and center 
sites on the Cu(lll) surface are shown in Fig. [2] their 
closeness illustrating the lack of corrugation on this sur- 
face, as in experimental findings [3] 0]. Similar results 
for bridge, atop, and center sites on Cu(100) show the 
PECs on the two surfaces to be very close to each other, 
just like for the experimental result [Hill]. The vdW-DF2 
equilibrium separation is about 3.5 A, like the value that 
is deduced from the experimental data and for reasonable 
physical assumptions about the parameters [3] |1] ■ 

A further refined comparison is provided by the bound- 
state eigenvalues for the point-of-gravity motion of H2 
on Cu(lll), both experimental data and those from the 
vdW-DF2 potential (Table |l|). In addition to PEC shapes 
and eigenenergy values, there should be comparison with 
values for well depth and equilibrium separation, exper- 
imental ones being (29.5 meV; 3.5 A) and those of the 
vdW-DF2 potential in Fig. [I] (37 meV; 3.5 A). Viewing 
the facts that (i) vdW-DF2 is a first-principles method, 
where characteristic electron energies are typically in the 
eV range, and (ii) the test system and results are very de- 
manding, as other popular methods deviate significantly 
more from the experimental curve (for instance, appli- 
cation of the DFT-D3(PBE) method gS], with atom- 
pairwise specific dispersion coefficients and cutoff radii 
computed from first principles, gives (— 88 meV; 2.8 A) 
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for the PEC minimum point). We judge this as very 
promising. So is the relative closeness of experimental 
and calculated eigenenergy values in Table [TJ 

The discrepancies between the eigenvalues signal that 
the vdW-DF2 PEC might not have the right shape for H2 
on Cu(lll). Reference ITU1 shows that vdW-DF2 bench- 
marks very well against the S22 data sets. It is pos- 
sible that the metallic nature of the H2/Cu(lll) sys- 
tem causes modifications in the details of the electrody- 
namical response; H2/Cu(lll) physisorption constitutes 
a very strong challenge for the density functional. For 
a well established conclusion, a more accurate theory is 
called for. 

In summary, accurate and extensive experimental data 
for the physisorption PEC of H2 on Cu(lll) are used 
to evaluate vdW-DF functionals and their adequacy for 
metal surfaces. The Cu(lll) surface is chosen here, as its 
flatness gives clarity in the analysis and eliminates several 
side-issues that could have made interpretations fuzzier. 
More generally, there exists an accurate data bank of 
experimental physisorption information that challenges 
every density functional approximation to produce rele- 
vant PECs. We propose that such surface-related PEC 
benchmarking should find a broader usage, supplement- 
ing, for example, S22 comparisons as an accelerator in 
density-functional development. 

Several qualitative similarities are found for both vdW- 
DF and vdW-DF2 functionals. The vdW-DF2 functional 
gives PECs in a useful qualitative and quantitative agree- 
ment with the experimental PEC, i.e. with respect to 
well depth, equilibrium separation, and curvature of PEC 
near the well bottom, and thus zero-point vibration fre- 
quency. This is very promising for applications of this 
nonlocal correlation functional at short and intermediate 
separations, relevant for the adsorption. However, the ac- 
curacy of experimental data is high enough to stimulate 
a more detailed analysis of all aspects of the theoreti- 
cal description. This should be valuable for the further 
XC-functional development. 
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